First order non-equilibrium phase transition and bistability of an electron gas 
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We study the carrier concentration bistabilities that occur to a highly photo-excited electron 
gas. The kinetics of this non-equilibrium electron gas is given by a set of nonlinear rate equations. 
For low temperatures and cw photo-excitation we show that they have three steady state solutions 
when the photo-excitation energy is in a certain interval which depends on the electron-electron 
interaction. Two of them are stable and the other is unstable. We also find the hysteresis region in 
terms of which these bistabilities are expressed. A diffusion model is constructed which allows the 
coexistence of two homogeneous spatially separated phases in the non-equilibrium electron gas. The 
order parameter is the difference of the electron population in the bottom of the conduction band 
of these two steady stable states. By defining a generalized free potential we obtain the Maxwell 
construction that determines the order parameter. This order parameter goes to zero when we 
approach to the critical curve. Hence, this phase transition is a non-equilibrium first order phase 
transition. 

Keywords: non-equilibrium phase transition; bistability; non-linear rate equations; semiconductor 
electron gas. 

Estudiamos las biestabilidades en la concentration de portadores que le ocurren a un gas de 
electrones altamente fotoexcitado. La cinetica de este gas de electrones fuera de equilibrio esta dada 
por un conjunto de ecuaciones de razon no-lineales. Para bajas temperaturas y fotoexcitacion en 
modo continuo (cw) mostramos que estas ecuaciones tienen tres soluciones de estado estacionario 
cuando la energfa de fotoexcitacion esta en un cierto intervalo, el cual depende de la interaction 
electron-electron. Dos de ellas son estables y la otra es inestable. Tambien, encontramos la region 
de histeresis en terminos de la cual estas biestabilidades son expresadas. Construimos un modelo 
de difusion que permite la coexistencia de dos fases homogeneas espacialmente separadas del gas 
de electrones fuera de equilibrio. El parametro de orden es la diferencia de la poblacion electronica 
en el fondo de la banda de conduction de estos dos estados estacionarios estables. Definiendo 
un potential libre generalizado del sistema obtenemos la construction de Maxwell que determina 
entonces al parametro de orden. El parametro de orden va a cero cuando nos aproximamos a la 
curva cn'tica. Por eso, esta transition de fase es una transition de fase fuera de equilibrio de primer 
orden. 

Descriptor es: transition de fase fuera de equilibrio; biestabilidad; ecuaciones de razon no-lineales; 
gas de electrones en semiconductores. 

PACS numbers: 72.20.Jv, 72.20.Dp, 78.20.Bh 



I. INTRODUCTION 

Bistability, threshold switching transitions, spatial 
pattern formation, self-sustained oscillations and chaos 
in semiconductors are related to non-equilibrium phase 
transitions P, These non-equilibrium phase tran- 

sitions in semiconductors have been studied in the 
past decades mainly in connection with the nonlinear 
generation-recombination mechanism, including impact 
ionization 0, d, HI- Optical and transport properties of 
semiconductors strongly depend on the electron popu- 
lation in the bottom of the conduction band and many 
experimental and theoretical studies on hot photo excited 
electron systems have shown that this electron popula- 
tion depends on the excitation energy [1,(1, 0). In this pa- 
per we present a study of the carrier concentration bista- 
bility that occur to an electron gas generated by a laser 
excitation and a study of a first order non-equilibrium 
phase transition between two homogeneous stable steady 



states of the excited electron gas. They both appear 
when the energy of the pump is varied. The order param- 
eter is the difference between the electron populations in 
the bottom of the conduction band of the stable states. 
It depends upon the effectiveness of the electron-electron 
interaction. We obtain the order parameter by defining 
a potential which allows us to make a construction simi- 
lar to the Maxwell construction for the equilibrium phase 
transition of a van der Waals gas[l, [§]. 

In following section we present the theoretical model 
we use to study a non-equilibrium electron gas in semi- 
conductors. In section third we deduce the stationary 
state solutions to the non-linear rate equations of the 
model and we make the estability analysis. By consider- 
ing a standard first order phase transition of the classical 
van der Waals gas, we analyze the non-equilibrium phase 
transition of the electron gas and obtain the correspond- 
ing Maxwell construction. We draw our conclusion in the 
final section. 
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II. THEORETICAL MODEL OF THE 
NON-EQUILIBRIUM ELECTRON GAS 

We describe the non-equilibrium electron gas in semi- 
conductors using the model equations given in Ref. 
Here, we give a brief description of these model equations 
and refer the reader to the details of their derivation to 
Refs. d, [H. Electrons in the bottom of the conduction 
band of a semiconductor play an important role in the dy- 
namics of the whole conduction electron gas. In general, 
electrons with an energy in excess less than the longitudi- 
nal optical (LO) phonon energy can not make transitions 
by emitting LO phonons. In the case in which the emis- 
sion of LO phonons is one of the dominant mechanisms, 
the nonequilibrium kinetics of the electron gas in the con- 
duction band of a semiconductor is given as follows [Toj. 
We define a set of energy levels, each one of them repre- 
senting an energy interval of width Ae of the conduction 
band. Although not strictly necessary, for simplicity Ae 
is set equal to fkuLO, the LO phonon energy 10]. We 
set the electron population \i m these energy levels and, 
based on the main interaction mechanisms, the nonlinear 
rate equations that give the temporal behavior of these 
populations are obtained. 

When an electron in level i emits a LO phonon it losses 
an amount of energy ftwLO and passes to the level i — 1 . 
If the frequency of this event is v£ then we have that 
the rate of change of population %i IS v t{Xi+t — Xi)- 
Similarly for the absorpion of a LO phonon we have the 
rate u~(xi-i ~ Xi) where v~ is the frequency of this 
process. The electron-electron interaction gives to the 
rate of change of the population Xi at level i two terms, 

ZN max Xtot(Xi+i - 2 Xi + Xi-i) and ZN max xo(Xi ~ Xi-i)- 
Their form come from considering the contribution to 
the rate of change of the population Xi of the interac- 
tion between electron populations of all the energy lev- 
els, and from the use of energy conservation [lOj. When 
an electron suffers a recombination it returns to the va- 
lence band and this gives us the term —wxi with w the 
frequency of the process. 

Then, by collecting all these contributions, we have the 
following set of rate equations which describes the kinet- 
ics of a photo-excited electron gas in semiconductors [lOj . 



dXi 
dt 



"oiXi+i - Xi) + v {Xi-i - Xi) 

+ZN max Xtot{Xi+\ - 2 Xi + Xi-i) 
+ZN max xo{Xi ~ Xi-i) + 9p$i,i P 



wXi 



(1) 



for i 7^ 0. For i = 0, since the emission of LO phonons 
by an electron is not possible, we have 



dxo 
dt 



V tX\ -! / „X0 + ^A r maxXtot(Xl - Xo) 

+ZN max xoXo + g P do,i p ~ wxo ■ (2) 



The electron populations Xi have been normalized to 
the maximum reachable electron concentration iV max and 
Xtot = X^Xi- The terms with the Kronecker delta is 



the generation contribution with generation rate g p . The 
main interaction mechanisms, the generation and recom- 
bination terms depend on the lattice temperature, carrier 
concentration and material parameters. 

Due to the electron-electron interaction, this set of rate 
equations is non linear. We should say that Eqs. (1) and 
(2) which are the main equations of our model came from 
a more general theory published in Ref. [l(|. See, also 
Refs. @,M EH for other details of the model equations 
and other applications. 

We make the following definitions v = v+ / 'ZN max , fi = 
is-/ZN max , u> = w/ZN max , x P = 9 P /ZN max , and e = 
e/Ae. By this normalization, we have eliminated one 
parameter. For low temperatures [lfl] . fj. <C v. Also, in 
steady state, Xv = w i Xtot = 1 and the total carrier 



concentration N = N n 



Hence, in this case, we have 



only two relevant control parameters, v and e p . Here 
e p = e p /Ae, the excitation energy in units of Ae. 

With these definitions and assumptions and consider- 
ing the finite differences as derivatives: 



Xi+i ~ Xt 
Xi+i - 2 Xi + Xi-i 

Xi Xi — 1 



dx 



de 
f_X 
de 



Ac 



dx 
de' 



Ae 



fx 

de 2 



(Ae) 



we have that the equations (TTJ) and that give the 
steady state of a cw photo-excited electron gas in semi- 
conductors are translated to the following differential 
equation 



d 2 x 



dx 



° = (1 ~ Xo) 1$ + {V + X0) ^ + ^ ~ £p) ~ ^ (3) 

The electron distribution function x(e, r) is continuous 
and positive definite in the whole interval < e < oo 
and has a discontinuity in its first derivative at e = e p . 
In addition, it must satisfy the conditions 



Xtot 



Xo 



de x(e, t) 



de x(e,r) 



(4) 
(5) 



The Eq. has the steady state solution^ 

s , , f Ae ae + Be-P £ ; < e < e„ 



where 



[{y + xS) 2 + Ml -Xo)] V2 T(^ + Xg) 
2(1 -Xt) 



(6) 



The electron population at the lowest level Xo 18 deter- 
mined self-consistently by the equation 



s W -ae e 



a + /3 1-Xo 



(7) 
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The coefficients A, B, and C are obtained using the con- 
tinuity in x(e), the discontinuity in its first derivative 
at e = e p , and Eq. (jU). In particular, we obtain that 
A/B = a//3 (See Ref. [7J for more details). 

To conclude this section we give the explicit expres- 
sion for the control parameter v that it is needed for 
the following discussion and we refer the reader to Refs. 
@, 0, El for more details and the expressions of other 
interaction mechanisms. Then, for v we have@ 



2^LO 



f 



(N q + 1) SujS^A (8) 



where T e is the effective electron temperature, fee is the 
Boltzmann constant, £ x and £ s are the static and optical 
dielectric constants, respectively, and ft is the Planck con- 
stant. Phonon population effects may be taken into ac- 
count in N q which is the phonon population at wavevec- 
tor of magnitude q. The screening in the electron-LO 
phonon interaction is taken into account in the factor [l3j 



( N 



where 



g 00 m(^ LO ) 3 
3 3 / 2 8ire 2 h 2 k B T c 



is the threshold value for the concentration in the conduc- 
tion band at which the screening becomes important (l3|. 
The electron effective mass and charge are m and e re- 
spectively. The screening in the electron-electron is given 
by the factor (loj 



in an easy way the relative importance of the electron- 
electron interaction against the electron-LO phonon in- 
teraction in terms of the carrier concentration and elec- 
tronic temperature. 

The electron-electron interaction which gives the non- 
linear character of the rate equation needs some clarifi- 
cation. We follow the ideas of Takenaka et aZ.[l4[ and 
Collet et al. [15j and use static RPA (Random Phase Ap- 
proximation) to obtain[9|, LLQ| 



ZN r , 



1 \fwmk~oT c 

^ 2 ft 2 £oo 



(9) 



The use of static RPA is justified for experiments which 
take place on longer time scales 16]. In particular, for 
bulk GaAs, under low temperature cw photoexcitation 
which creates a not very high carrier concentrations (less 
than 1 x 10 18 cm -3 , where many body and occupation 
effects are not important and the electron-LO phonon 
and electron-electron dominates the dynamics of the car- 
riers) the equations of the model are expected to be 
valid [HI E3]- Degeneracy effects are also negligible for 
concentrations below 1 x 10 18 cm" 3 [17]. As we men- 
tion before the equations of the model came from a 
more general theory published in Ref. [13]. These gen- 
eral equations take into account degeneracy effects and 
can be generalized to take into account quantum ef- 
fects, like exch ang e, by using more appropiate scatter- 
ing frequencies [18j . However, quantum effects are ex- 
pected to be important in very short times and very small 
distances 



be imp 



See = 1 



N 

and becomes important when the carrier concentration 



III. STATIONARY STATES AND STABILITY 
ANALYSIS 



A PO = 



4771500 (fc B T c 
ir 2 ft 2 e 2 



The frequencies associated with the collision mecha- 
nisms are in general dependent upon the energy. There- 
fore, v is a function of e but a smooth oneQ . The analysis 
can be considerably simplified if we average it over the 
conduction band. The factor A given by 



A = 



£m ax i 

ds—= In 



We' 



1 



1 



is this averaging process and e max is the maximum en- 
ergy in the conduction band that can be reached by an 
electron [To| . 

The expression for v, Eq. ©, is one of the simplest ap- 
proximations we can make which allows us to determine 



In Fig. 1 we show the low temperature dependence 
of Xo as a function of the control parameters e p and v. 
Fig. la shows the dependence of Xo upon the control 
parameter v for various values of e p . We have taken 

uj = 0.02 for the recombination. When e p ~ 16 and 
the efficiency of the electron-electron is increased (y is 
diminished) we found a bifurcation, a range of values of 
v where we have more than two solutions for Xo- 

The dependence of \ s upon the control parameter e p is 
shown in Fig. lb for several values of v. Notice that, for 
example, the curve with v — 0.02, for e p < e p < e p 2 ^ has 
three possible values of Xo, that are labelled Xo \ Xo j 

(3) 

and Xo ■ We will show that two of them are stable and 
the other one is unstable. In addition, for v c w 0.2275 
we have the critical curve which separates curves which 
has regions of e p with three possible solutions for x§ from 
curves with just one solution. The critical point is given 
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FIG. 1: Xo as a function of the control parameters v (a) and 
e p (b). The number on each curve is the value of e v (a), v 
(b). 



by £p, Xo an d v c that satisfy 



de r . 



dXo 
\dx 2 o 



= 
= 



(10) 

(11) 



We must point out that v = z/+ /ZN max , then Fig. lb 
can be seen as a family of curves in which the electron- 
electron interaction is changed, for example, by changing 
the electron concentration. For GaAs bulk semiconduc- 
tor, at 4 K, using N = 4.1 x 10 17 cm" 3 and T = 300 
K in Eq. ([8]) we obtain v = 0.02. The critical curve is 
obtained for N = N c « 6.11 x 10 16 cm" 3 . 

The extrema of the curve £ P (xo) f° r a g lven v are e P 

( 2) 

and e p respectively. They are found by solving the equa- 
tion (d£ p /dxo)v = 0- They are the boundaries of what is 
call the hysteresis region. If we increase e p from a value 

less than Sp 1 ^, we will move until the maximum value e p 2 ^ 
is reached. Then an abrupt jump will take place from a 

(2) 

Xo value at the maximum e p to the value marked with 
the arrow at the left of Fig. lb. As we decrease e p from 

(2) 

a value greater than e p we will move until we reach the 



FIG. 2: Hysteresis region (hatched) in the phase space (v,e p 
for two values of lu. 



value of Xo at the minimum. Then the system will un- 
dergo an abrupt jump to the value of Xo marked with the 
right arrow in Fig. lb. In Fig. 2 we show this hysteresis 
region for two values of the recombination parameter ui. 

Let us consider the perturbation to the steady state, 
x(e, t) = x s ( £ ) + &x{ e i T ) produced by a time dependent 
excitation, Xp(t) — X P + ^XpWj where x p is the constant 
excitation that produces the steady state. Notice that 
Xtot and xo ar e time dependent and they are given by 
the relations 



Xtot = 1 + Sxtot 
Xo = Xo + s Xo 



(12) 



From Eq. ((5J) for xo an d assuming that e p > 1, we have 
for its steady state 



VX\ + Xt -Xo + Xo - ^Xo = 



(13) 



and the time evolution of its perturbation Sxo is given 
by 



— 8xo = v5xi + (Sxi ~ *Xo)(l +<^Xtot) +<^Xtot(Xi - Xo) 



+2xo<5xo + (<5xo) 2 - wJxo 



(14) 
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which reduces to first order 



-7-^X0 = v8xi + Sxi ~ <^Xo + 2xo<5xo - u5\a 
dr 



S Xi ,8X1 , , s 
v- Y 1 + 2x - w 

<>Xo £>Xo 



<5xo (15) 



We now define the function ip(xo) that allows us to 
perform the stability analysis. From Eq. (fT3|) 



■0(Xo) = vx\ +Xi - Xo +Xo - w Xo (16) 
and the steady state condition (Eq. (fT3])1 becomes 

^(XS) = (17) 
The derivative of V'(Xo) with respect of xo is 

<¥(xo) dxi d X i , , 

— = v- V- l + 2xo-w (18) 

dxo dxo dxo 

Now, in steady state xi is function of Xo an d the pertur- 
bation 5xi is written as 



Therefore, 



I dxi 

\dxo 



d ^ d^? 
dr dxo 



and the solution for <5xo is 



5xo 



Xo=Xo 



5xo 



(19) 



Xo=X5 



^Xo(t) = 5xo(0)exp 



A steady state, given by a Xo> is stable if 



" dtp 




dxo 


T 

Xo=Xg 



X0=XS 

Eq. ([7]) can be written as 



< 



(20) 



(21) 



1 



In 



(1 - Xo)x s o 



a a + j3 



-0 



(22) 



In steady state Eqs. (|16|) and (|22|) are equivalent, there- 
fore 



<Kxo) 



-iln 

a 



a a + f3 
(!-Xo)Xo-— — 



(23) 



In Fig. lb we plotted function of Xo which means 

that we plotted the function ip(xo) and this establishes 
the stability of the states Xo~\ Xo 1 an d xif' '• The middle 
root is an unstable steady state, while Xo an d Xn are 
stable steady states. 



IV. FIRST ORDER NONEQUILIBRIUM PHASE 
TRANSITION AND MAXWELL 
CONSTRUCTION 

A thermodynamic system remains homogeneous and 
stable if the criteria of instrinsic stability is satisfied [l!|, 



dP 
dV 



^7 <0 



(24) 



where P is the pressure of the system, V is its volume 
and T is its temperature. When this condition is violated 
a phase transition ocurrs(l9j. 

Then, from Eq. (|21j) , we have the following correspon- 
dence [23 



Xo 

Sp 

v -> T 



V 

p 



We see that Eq. ( 122)) turns out to be the equation of 

state. The homogeneous states Xo anc i Xo^ are nonequi- 
librium stable steady states of the system. Therefore, we 

call the transition between states Xq anc i Xq 3 \ a nrs t 
order out of equilibrium phase transition. The order pa- 
rameter is Xo 3 "* — Xo^ an d still remains unknown. 

So far, we have assumed that the populations xo> Xi> 
and so on, are constants and homogeneous in space. Now, 
we assume that there exist spatial inhomogeneities, xo = 
Xo(r), which produce spatial gradients, and diffusion of 
this population. We also assume, from the structure of 
Eq. ([3]), that the spatial and temporal behavior of the 
populations Xii i > 0, is given through xo- Then, we 
have the following equation 



dxo 
dr 



fxo 
dz 2 



(25) 



We suppose, for simplicity, that xo depends only on the 
spatial coordinate z. Here, k is the xo difussion constant. 

Let us introduce the "potential" $(xo) with the 
definitionji, i] 



d 

V'(xo) = -tt—Hxo) 
oxo 

Then, the steady state satisfies the equation 
d 2 Xo 



dz 2 



■ o — ^(Xo) 
dxo 



(26) 



(27) 



We have seen from Fig. lb that, when v — 0.02 and e p 

is in the energy interval, Ep < e p < e p 2 \ the system 
has two homogeneous stable steady states. Let us find a 
solution Xo{ z ) such that, xo(+°°) = Xo an d Xo(~°°) = 
xb 3%> ■ In such a case two steady states coexist. Obviously, 
the "potential" < I > (xo) has two maxima in Xo an d X^ ■ 
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Coexistence ocurrs for a value of e p such that the two 
maxima are indistinguishable for the system, 

2(X$ ) ) = *(X$ ) ) (28) 

Then 

(3) 

= $( x ^) - = [ ° dxo^(xo) 

J Xo 

= / dxo (vXi + Xi - Xo + Xo - wxo) (29) 

The last equation is the Maxwell construction for the va- 
por pressure in the van der Waals gas from which the 
order parameter Xo — Xo can be calculated. We also 
notice that — $ corresponds to the Hemholtz free poten- 
tial. Moreover, the condition that $ is at a maximum 
in a stable steady state corresponds to that the gener- 
alized free potential H — — $ is at a minimum. This is 
consistent with the condition that an equilibrium ther- 
modynamic _system is in a state of minimum Hemholtz 
free energy[19l]. 



study them in terms of the hysteresis region. The dif- 
fusion model consider here allows the coexistence of two 
phases of the open far from equilibrium gas. The phases 
are distinguished by the two values of carrier concentra- 
tion at the bottom of the conduction band, namely, Xo^ 

(3) 

and Xo ■ These two corresponding homogeneous stable 
steady states are obtained using the Maxwell construc- 
tion. We introduced the function ip and the stability 
condition for these states was given by Eq. (|2"Tj). A gen- 
eralized free potential was defined by Eq. (|26|) from which 
we obtained Eq. (j2"5)) that corresponds to the vapor pres- 
sure Maxwell construction of a van der Waals gas. The 
difference of the electron population in the bottom of the 
conduction band is the order parameter and is then cal- 
culated and goes to zero when v approaches the critical 
value v c ss 0.2275. This is the reason we call this phase 
transition a first order nonequilibrium phase transition. 
For bulk GaAs at 4 K the critical curve is obtained us- 
ing N = N c fa 6.11 x 10 16 cm" 3 and T c = 300 K in 
Eq. {§]). We have to point out that the electron-electron 
interaction, which gives the nonlinear character of the 
rate equations, is the necessary main ingredient for the 
existence of the phase transition. 



V. CONCLUSIONS 



We have found carrier concentration bistabilities in 
a low temperature nonequilibrium electron gas and we 
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